Geometry-related magnetic interference patterns in long SNS Josephson junctions 
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We have measured the critical current dependence on the magnetic flux of two long SNS junctions 
differing by the normal wire geometry. The samples are made by a Au wire connected to W contacts, 
via Focused Ion Beam assisted deposition. We could tune the magnetic pattern from the monotonic 
gaussian-like decay of a quasi ID normal wire to the Fraunhofer-like pattern of a square normal wire. 
We explain the monotonic limit with a semiclassical ID model, and we fit both field dependences 
with numerical simulations of the 2D Usadel equation. Furthermore, we observe both integer and 
fractional Shapiro steps. The magnetic flux dependence of the integer steps reproduces as expected 
that of the critical current /c, while fractional steps decay slower with the flux than Ic- 



PACS numbers: 74.45. +c, 73.23.-b, 74.40. Gh 

Introduction - A non-dissipative supercurrent can be 
transmitted between two superconductors (S) through a 
nanometer-thin insulating layer (I), when a phase differ- 
ence is imposed. Superconducting correlations can also 
penetrate into a micrometer-long, non-superconducting 
coherent metal (N) over lengths greater than ^5, the su- 
perconducting coherence length in N. A supercurrent can 
then flow through long (L > ^5) SNS junctions, provided 
that the phase coherence is preserved in the normal wire. 
The supercurrent amplitude only depends on the normal 
metal length L and diffusion coefficient D (which set the 
characteristic Thouless energy Eth = hD/L'^) and on 
the normal metal resistance.^ The supercurrent thus re- 
flects the transport mechanisms in the normal wire, and 
is affected by the interference and diffraction phenomena 
present in the normal metal as a result of the phase co- 
herence. 

We have studied the behavior of the supercurrent in long 
SNS junctions in a perpendicular magnetic field. In a 
previous work,!^ we observed that the maximum super- 
current, the critical current /c, monotonously decreased 
with the magnetic field following a quasi-gaussian depen- 
dence. This behavior is different from the interference 
Fraunhofer pattern usually observed, for example, in SIS 
junctions, short SNS junctions (L < ^5), magnetic SFS 
Josephson junctions, etc. This difference is induced by 
the different aspect ratios of short weak links and long 
SNS junctions. Indeed, in SIS junctions the thickness 
of the junction is limited to a few Angstroms, to permit 
the tunneling of Cooper pairs; similarly, in SFS junctions 
the magnetic layer thickness has to be shorter than a few 
tens of nanometers, so that Cooper pairs are not broken 
by the internal exchange field. In contrast to the wide 
and short SIS and SFS junctions, SNS junctions offer the 
interesting possibility to explore a broad range of aspect 
ratios, since the length of the normal metal is only lim- 
ited by the phase coherence length, which can be as long 



as a few microns at low temperature. In this paper we 
explore different geometries of long SNS junctions. Both 
monotonic and non- monotonic Ic{B) dependences have 
been observed before, but we show for the first time that 
we can tune the Ic{B) curve from an interference pattern 
to a quasi-gaussian monotonic dependence by varying the 
normal metal's aspect ratio. 

Monotonic Ic{B) dependences have been observed in 
ballistic long SNS junctions, formed by a normal bidi- 
mensional InAs electron gas connected to supercon- 
ducting Nb contacts larger than the London screening 
length.^ Their magnetic field dependence resulted from 
the screening currents in the Nb. On the contrary, the 
diffusive junctions investigated in the present work are 
contacted by thin disordered superconducting wires in 
which the magnetic field screening is negligible. We show 
that the geometry dependent magnetic field decay can be 
explained taking into account only the interferences be- 
tween Andreev pairs' trajectories in the normal metal. 
As a reminder, we first consider the case of a wide short 
SIS junction. A magnetic field in the SIS junction plane 
penetrates in the insulating layer of thickness d and in 
the superconductors nearby over a length A^, the Lon- 
don penetration length. In a magnetic field B = —Bz of 
vector potential A = Byx, the phase shift of the Cooper 
pairs tunneling at different points of the junction width 
is (Fig. [Ijc)): 

e{y) = 27r^ / A,dx = 2^-^ (1) 

where ^{y) is the flux through the surface Sy = {d -\- 
2Al)^, and ^0 = h/{2e) is the quantum flux. The cur- 
rent is obtained by integrating over the junction surface 
the supercurrent density j = jcSin{S + ^), taking into 
account both the superconducting phase difference be- 
tween the contacts S and the phase due to the vector po- 
tential. The critical current dependence on the magnetic 
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flux ^ = B{d-\-2XL)w, follows the well-known Fraunhofer 
pattern, a diffraction pattern created by the interference 
between the ballistic trajectories over the junction width: 



Ic = /c(0) 



^0 

7r$ 



sm 



^0 



(2) 



If we now consider SNS junctions, we expect strong 
differences between junctions containing a short and 
wide or a long and narrow normal wire: in a short 
wide wire the phase difference between the trajectories 
comes from the phase distribution along the junc- 
tion width, while in a narrow long wire, the phase of 
each trajectory is accumulated along the junction length. 

The samples - We have fabricated long SNS junctions 
where a Au normal wire links two superconducting W 
contacts. First, a 50 nm thick Au wire is drawn by e- 
beam lithography and deposited onto a Si02 substrate. 
We use 99.9999% pure gold, with a content in magnetic 
impurities (Fe) smaller than 0.1 ppm. This insures a 
long phase coherence length ~ 10 fim below 50 mK, 
measured in a separated weak localisation experiment. 
We then deposit the superconducting contacts in a Fo- 
cused Ion Beam (FIB). After slightly etching the Au wire 
with the FIB to remove possible impurities on its sur- 
face, we inject a metallo-organic vapor of tungsten hex- 
acarbonyl over the sample. This vapor is decomposed 
by the focused Ga+ ion beam, and a disordered W al- 
loy is deposited on the substrate. The wires produced 
are composed of tungsten, carbon and gallium in vary- 
ing proportions (in our case, the atomic concentrations 
are roughly 30% W, 50% C and 20% Ga)."^ The super- 
conducting critical temperature of the wires produced is 
Tc ~ 4K, an order of magnitude higher than the bulk Tc 
of W."^ This could be due to the inclusion of Ga, which is 
itself a superconductor with Tc = IK. The critical mag- 
netic fleld of the wires is also strikingly high: at IK, 
He = 7T.^ The W wires are 200 nm wide and 100 nm 
thick. The dependence of the superconducting properties 
of these wires on the deposition conditions have been in- 
vestigated in detail in W. Li et al.J The superconducting 
gap as well as the Abrikosov flux lattice have been studied 
by STM experiments.!^ The investigation of proximity in- 
duced superconductivity in metallic nanowires contacted 
by FIB has also been recently performed.^ The advan- 
tages of this technique are the deposition of the material 
of practically any shape and size, without any mask, and 
the good quality of the interface created. The main disad- 
vantage is the Ga contamination of about 250 nm around 
the deposited wires. The long SNS junctions created by 
FIB-assisted deposition are comparable to junctions cre- 
ated with more standard fabrication methods. We thus 
recovered the general results for the voltage vs. current 
curves and the temperature dependence of the critical 
current .^13 

To investigate the influence of the geometry on the /c(^) 
dependence, we have designed two samples with different 
aspect ratios: sample WAu-Sq is 1.2 jam long and 1.75 
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FIG. 1: (color online) Scanning Electron Microscope images 
of samples WAu-Sq (a) and WAu-N (b). One can see as a 
clear halo the W contamination around the superconducting 
contacts, (c) Scheme of a SIS Josephson junction, (d) Scheme 
of the ID model developped. 




FIG. 2: (color online) Sample WAu-N normalised critical cur- 
rent vs. normalised flux at T=60 mK. Red line: Gaussian 
function with cr=0.5. Yellow line: numerical simulation of 
the 2D Usadel equation, where the flux has been rescaled by 
a factor 2.5. Inset: comparison between a Gaussian depen- 
dence (dotted line), the semiclassical model prediction for a 
ID normal wire (black line), the analytical result of Usadel 
equation in the limit L » w (red line) and the numerical 
simulation of the 2D Usadel equation for the aspect ratio of 
sample WAu-N (blue line). /c,maa;=2.4 /xA. 



Ilia wide, with an aspect ratio L/w = 0.7, while sample 
WAu-N is 1.53/im long and 0.34 /im wide, with an aspect 
ratio L/w = 4.5 (Fig. [lja),(b)). 

The /c(^) normalised curves for samples WAu-N and 
WAu-Sq are shown in Fig. |2] and |3] respectively. They 
illustrate the important role of the aspect ratio: sample 
WAu-N displays a quasi-gaussian decay of the critical 
current, while sample WAu-Sq displays oscillations 
which recall a Fraunhofer pattern. 

Semiclassical model - To explain the behavior of sam- 
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FIG. 3: (color online) Sample WAu-Sq normalised critical 
current vs. normalised flux (blue dots). A trapped flux of 
4.3 G has been subtracted. Red line: numerical simulation 
of the 2D Usadel equation for a junction with aspect ratio 
L/w—0.7 and W wires inductance £=11.5 pH. Inset: raw 
data for Ic{B), with Ic,max =38.8 /iA. T=60 mK. 



pie WAu-N, a semi-classical model was developed.™ The 
aim is to model a long diffusive SNS junction, with 
w « L << L(j)^LT (ID geometry), in a perpendicular 
magnetic field B = —Bz of vector potential A = Byx. 
We suppose a sinusoidal current-phase relation for sim- 
plicity: j = jc sm{5-\-0ij). Oij is the phase acquired along 
the trajectory starting at point i and finishing at point 
j. S is the phase difference between the two supercon- 
ducting contacts, and is independent of the position of i 
and j, since in this long ID geometry one can ignore the 
dephasing along the junction width in comparison to the 
one along the junction length. The total current is the 
average over all possible trajectories Cij: 
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Considering a reference trajectory 1-2 with a phase dif- 
ference we have Oij = Oa + AOij (see Fig. [^d)). The 
critical current is the maximum of the current /: 
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Given the great number of possible trajectories, the cen- 
tral limit theorem sets the distribution of the trajectories 
length to a Gaussian distribution. The phase asso- 
ciated to the trajectory zj, follows then also a Gaussian 
distribution. Thus: 



-l/2<(A^,,)2>c^^. 
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the surface Sij defined by the trajectory ij (Fig. [TJd)): 
2e 



A.j^ dx / A.j^ dx 



- f A^dx = 2n— 
(6) 

Introducing eqj6]in eqjS) we find that the critical current 
decays as a Gaussian function: 
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where = [S^ / < Sfj >c,J/(47r2), and 5' = ^ x L. 
If the trajectories are ballistic, a Gaussian decay with 
(7=0.28 is found.[I2 The Gaussian fit of sample WAu-N 
gives (7=0.5, while in our previous results,^ we obtained 
(7 in the range 0.74-1.8 for aspect ratios L/w in the range 
1.9-10.4. 

The exact calculation describing diffusive trajectories 
from point i to point j gives: 



Vs ^0 



sinh ( ^-t- 



(8) 



where AOij = Oij — Oa is proportional to the flux ^ij in 



The difference between the above distribution and a 
Gaussian one is shown in the inset of Fig. [2] While this 
model explains the origin of the quasi-gaussian shape of 
our /c(^) curves, the experimental decay is faster than 
the predicted one. A more precise calculation, taking into 
account the 2D nature of both samples, is thus necessary. 

2D Usadel equation - J.C. Cuevas and F.S. Bergeret!^ 
have solved the 2D Usadel equation for long diffusive 
SNS junctions with low-resistance interfaces and for 
different aspect ratios of the normal wire. A complete 
field penetration in the normal metal {w < Xl,n) is 
assumed, neglecting any Josephson current screening 
effects. Moreover no inelastic scattering is considered. 
The authors found two very different limits: 

- w « L : in this ID limit, the field acts as a pair- 
breaking mechanism. The critical current is then mono- 
tonically reduced: at T=0 Ic{B)/Ic{0) - g-^-^^s (<^>/<^>o)' 
(inset of Fig. [2|. 

- w >> L : in this case, the main effect of the field is 
to modulate the phase over the junction width, giving 
rise to an interference pattern, which is identical to 
a Fraunhofer pattern for small enough aspect ratios 
{L/w <0.04). 

The authors do not suppose a purely sinusoidal current- 
phase relation to start with, but calculate the complete 
I{S) relation of a long SNS junction. However, the wide 
junction limit corresponds exactly to the Fraunhofer 
pattern obtained in the case of a sinusoidal current-phase 
relation. 

Experimental results - In agreement with the predic- 
tion of J.C. Cuevas et al. for the aspect ratio I//k;=4.5, 
junction WAu-N does not exhibit any oscillating pattern 
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(Fig. |2|. The field decay cannot however be fitted by 
the numerical simulation for the corresponding flux scale 
$/^o- the experimental decay is faster than the pre- 
dicted one by a factor 2.5. This can be explained by 
taking into account a non-perfect interface, with an in- 
terface resistance Ri given roughly by Ri = 2Rn ~ 5^1 
(the ratio Ri/Rjsj = 2 leads indeed to a rescaling of the 
perfect interface curve by a factor 2.3). This non-perfect 
interface may be attributed to the right W contact, which 
is just at the extremity of the gold wire (see the SEM im- 
age 0b). 

To understand the experimental Ic{^) of sample WAu- 
Sq, it is necessary to take into account not only the as- 
pect ratio of the normal wire, but also the flux correc- 
tion due to the inductance of the W contacts, evident 
in the tilt of the central peak of the Ic{^) curve. The 
'internal' flux ^int through the normal wire is thus the 
sum of the applied external flux <l> and the flux created 
by the current flowing in the superconducting contacts, 
= C X Ici^int)^ following the below relation. 

^int=^^CIc{^int) (9) 

The experimental curve should then be com- 
pared not to the calculated curve PJ^{^int) but to 
jth^^^ = If{^^nt - (^,nt)). The bcst fit of our 
data for sample WAu-Sq according to the above expres- 
sion yields £=11.5 pH. This value is compatible with 
the kinetic inductance, much larger than the geometrical 
one, of similarly made W wires. 

When comparing our experimental results with the nu- 
merical simulation of 2D Usadel equation for L/w = 0.7 
modified by the self-inductance, we find a good agree- 
ment in the position of the zeros for <l> = 5 with 
S = 3.25/im^. This corresponds roughly to the whole 
surface of the normal metal square, slightly larger than 
the surface between the contacts. The amplitude of 
the measured oscillations following the central peak 
decreases slower than predicted, and the minima of the 
first periods do not go to zero; this can be qualitatively 
explained as the effect of a non uniform current distri- 
bution in the normal metal.*^ 

Shapiro steps - We have measured the differential 
resistance dV/dl when irradiating the junctions with 
microwaves from an antenna. We have observed Shapiro 
steps, in the form of microwave- induced zero resistance 
dips at V = ^nf. They result from the resonance of 
the ac voltage at frequency / induced by the microwaves, 
with the current oscillations at frequency fj = 2eV/h^ 
due to the ac Josephson effect at finite voltage V (Fig. 
|4]^a),(b)). In addition to the Shapiro steps at integer 
n, we find fractional Shapiro steps at n=l/2, n=l/3, 
n=3/2 and n=l/4. Both integer and fractional steps can 
be observed at voltages larger than the Thouless energy 
{Eth/^ = 5.3 /iV), and for frequencies larger than the 
minigap {A/h = 4 GHz); the voltage was directly de- 
duced from the measured V(I) curves. Fractional Shapiro 
steps appear in SNS junctions as a consequence of a 



non-sinusoidal current-phase relation. The additional 
harmonics in the current-phase relation can be generated 
by multiple coherent Andreev reflexions (MAR), when 
the coherence length is much longer than the N length, 
or by non-equilibrium effects. Fractional Shapiro steps 
reflect the behavior of each harmonic individually: the 
step at n=l/2, for example, is generated by the second 
harmonic and is proportional to its amplitude, having 
the same dependence on field and temperature. We have 
studied the magnetic field dependence of both integer 
and fractional steps amplitudes. As expected, the field 
dependence of the integer Shapiro steps follows roughly 
that of the zero-th order step, the critical current (Fig. 
[4jc)). This is not the case however for the fractional 
steps: in Fig. [4] we show that the decay with the 
normalised flux of step n=l/2 is slower than that of the 
critical current. If this fractional Shapiro step was due 
to the MAR at equilibrium, the field dependence would 
show a periodicity half that of the critical current, cor- 
responding to the double length covered by the Andreev 
pairs, and the zeros of the Fraunhofer pattern should be 
found at multiples of <^>o/2. This is evidently not the 
case, so that the supplementary harmonic in the current- 
phase relation may be traced back to non-equilibrium 
effects, as already suggested in P. Dubos et al)^ and F. 
Chiodi et a/.,^^ where the effect of the magnetic field at 




FIG. 4: (color online) (a) Differential resistance dV/dl of sam- 
ple WAu-Sq in presence of an irradiation at frequency f=2.8 
GHz. (b) Dc voltage vs. Shapiro step order for sample WAu- 
Sq; lines show the predicted dependence V — h/2enf. (c) 
Normalised flux dependence of Shapiro steps current ampli- 
tude for n=l and n=l/2 (symbols) compared to the critical 
current dependence under microwaves, which explain the de- 
crease of its amplitude (continuous lines); the trapped flux 
of 4.3 G, corresponding to 0.7 $o is partially compensated by 
the inductance-related tilt of Ic{B) central peak. 
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the lowest temperatures was shown to increase the to- 
tal critical current in the out of equilibrium SNS junction. 

Conclusions - We have measured the Ic{^) curves 
for two different geometries of long SNS junctions. 
The samples are made by a Au wire connected to W 
contacts, via FIB-assisted deposition. We have observed 
a monotonic gaussian-like decay for a quasi ID normal 
wire, in contrast to the Fraunhofer-like interference 
pattern of a square normal wire. We explain the 
monotonic limit with a semiclassical ID model, and fit 
both field dependences with numerical simulations of 



the 2D Usadel equation. Moreover, we have observed 
both integer and fractional Shapiro steps and their 
dependence in magnetic field. While integer steps 
follow as expected the field dependence of the critical 
current, fractional steps decay slower than /c- This is 
incompatible with equilibrium MAR-originated steps, 
but may be explained as an out-of-equilibrium effect. 
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